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^ ! Abstract 

^~~) ' In a sample of 9.66 x 10 BB pairs collected with the CLEO detector we 

■ make the first observation of B decays to an r]c and a kaon. We measure 

^ ! branching fractions B{B+ VcK+) = (0.69l[]:|? ± 0.08 ± 0.20) x 10"^ and 

> ; B{B° r]cK°) = {imtoil ± 0.12 ± 0.31) x lO'^, where the first error is 

I statistical, the second is systematic and the third is from the r/c branching 

' fraction uncertainty. From these we extract the r/c decay constant in the 

^ ■ factorization approximation, /^^ = 335 it 75 MeV. We also search for B decays 

Q . to a XcO and a kaon. No evidence for a signal is found and we set 90% CL upper 

O ; limits: BiB+ XcoK+) < 4.8 x lO""^ and B{B° XcoK°) < 5.0 x lO""^. 
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Two-body B decays to a charmonium state and a kaon have recently received substantial 
attention because of their importance for studies of CP violation in determining the angles 
of the CKM unitarity triangles and because of the observation of the unexpectedly large B 
decay rate to ri'X Among several theoretical explanations for the latter, a substantial 
intrinsic charm component of the rj' has been proposed 0,^. If this is the case, the t]' 
can be produced by the axial part of the b ccs{d) process, which also produces the rjc- 
Exclusive B decays to charmonium states are also of theoretical interest as a testing ground 
for the QCD calculations of quark dynamics and factorization. In the absence of enhancing 
mechanisms, the B decay rate to rjcX is expected to be comparable to that for the B decay 
to J/tpX Experimentally, little is known about rjc production in B decays. The only 

pubhshed result is from a 1995 CLEO study 0], which used 2.02 fb^^ of data collected at 
T(4S') and obtained an upper limit on inclusive rjc production: B{B Vc^) < 0.9% at 90% 
CL. 

The color-singlet production of Xco in B decays vanishes in the factorization approxima- 
tion as a consequence of spin-parity conservation. However, the color-octet mechanism P,p!0 



allows for the production of the Xco P-wave 0^+ state via the emission of a soft gluon. No 
information on B decays to Xco is available at present ||13|| . 

In this Letter we report results from the analysis of 9.13 fb~^ of e^e^ annihilation data 
collected with the CLEO detector at the Cornell Electron Storage Ring (CESR), taken at 
the T(4S') energy, corresponding to 9.66 x 10^ produced BB pairs. In addition, 4.35 fb~^ 
of integrated luminosity were taken 60 MeV below the T(4S') resonance in order to study 
backgrounds from light quark production (referred to as continuum). 

The data were taken with two configurations of the CLEO detector, called CLEO II and 
CLEO II. V. In the CLEO II configuration of the detector [|1^, charged particle tracking is 
provided by three cylindrical drift chambers immersed in axial solenoidal magnetic field of 1.5 
T. Charged particle identification (PID) is made possible by a time-of-flight system (TOP) 
outside of the outermost tracking chamber and by the measurement of specific ionization 
loss (dE/dX) in the tracking system. Photon and electron identification is provided by 
a high resolution electromagnetic Csl (Tl) calorimeter. The muon system is the outermost 
subdetector consisting of three superlayers of wire counters interspersed with steel at different 
absorption lengths. The CLEO II. V configuration differs from the CLEO II configuration in 
that the innermost drift chamber was replaced by a three layer double-sided silicon vertex 



detector |T5| and that a helium-propane gas mixture, instead of argon-ethane, was used in 
the main drift chamber. These changes led to improved momentum and dE/dX resolution. 

We reconstruct the rjc in the decay modes rjc ^ (fxp ^ K~ K~ and rjc K^K^tt^ 
Q. The XcQ is searched for in its decay modes Xco K^K^ and Xco ~^ vr+vr^. For calculation 
of efficiencies, we use the branching fractions i3(xco ^ K^K") = (0.586 ± 0.086)% and 
^iXco 7r+7r~) = (0.496 ± 0.066)% obtained by averaging the PDG values [0 with the 
recent BES results [|I6|. The detector simulation is based upon GEANT 3 p7|. 



Candidate primary tracks must be well measured and come from the event vertex. Neutral 
kaons are identified as a tt+tt^ pair coming from a displaced vertex. The mass resolution is 



Charge conjugation is implied throughout this Letter. 
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3.6 MeV and we select events within 8 MeV of the mass A phi meson candidate is 
selected as a K^K~ pair in the mass window 1.00 < M{K'^ K^) < 1.04 GeV. 

We reconstruct the B mesons by combining an rjc or a Xco with a charged or neu- 
tral kaon. The candidate events are identified by using the difference between the re- 
constructed and beam energies, AE = E{B) — Eheam, and the beam constrained mass, 



M{B) = ^^Lm - \Pb?- The resolution of AE is about 17 (15) MeV for CLEO II (CLEO 
II. V). The uncertainty in M{B) is about 2.6 MeV and is dominated by the beam energy 
spread. Events for which multiple combinations pass the selection criteria are assigned a 
weight equal to the inverse of the number of candidates passing the selection. The average 
number of candidates per event is about 1.3 for decays in the rjc KgK^i^^ submode and 
is less than 1.1 for all other channels. 

To minimize the combinatorial background coming from other B decays and from con- 
tinuum production we impose PID criteria. For the 7]c 00 submode we require that the 
charged kaons have dE/dX and TOF measurements within 3 standard deviations (a) of the 
expected values, when such measurements are present. The r/c — > K^K^n'^ submode has 
more background and the PID consistency requirements in this case are at a more stringent 
level of 2o", and at least one of the PID measurements has to be present for the charged 
kaon candidates. The secondary tracks from the Xco decay have momenta between 1.0 and 
2.7 GeV, where there is little PID separation between kaons and pions, therefore, no PID 
requirements are imposed. 

Most of the background comes from continuum production, which at CLEO has jet-like 
characteristics. We minimize this background by employing the ratio R2 of the second and 



zeroth Fox- Wolfram moments |jT8|. For isotropic events, R2 is nearly zero and for jet-like 
events it is close to one. We select events with R2 < 0.25 for decays with rjc and R2 < 0.3 
for decays with Xco- In addition, we impose a lepton veto on the bachelor kaon candidate to 
remove possible contamination due to B semileptonic decays. 

We extract the signal yield by using an unbinned extended maximum likelihood fit 
method described in Ref. |jT9|. For the B rj^K analysis, the fit variables are the beam 
constrained mass M{B), the energy difference AE, the rjc candidate mass and the cosine 
of the angle between the direction of the B candidate and the beam axis cos{6b)- For the 
rjc — > (f)(f) mode we also use the angle x between the planes formed by the kaons from the 
decays in the rjc rest frame . 

The signal M{B) and AE probability density functions (PDFs) are parameterized by 
a Gaussian and a sum of two Gaussians with the same mean, respectively. The rjc mass 
is represented by a Breit-Wigner function convolved with a sum of two Gaussians with the 
same mean. We use 13.2t3;2 MeV for the r]c width [0. The signal shape is expected to vary 



as sin {6b) and sin (x). The background shape is represented by an end-point function (the 



product of a 2nd degree polynomial in y 1 — {M{B)/Ey and a phase-space factor) in M{B) 
and is linear or constant in all other fit variables. The parameters of the PDF shapes are 
extracted from the Monte Carlo simulation of the signal and backgrounds. We combine the 
yields from different data sets and from different 7]^ submodes by adding the log-likelihood 
functions. The yields and efficiencies are given in Table |. 

Systematic uncertainties from the modeling of the PDF shapes are included in the fit 
result by varying the shape parameters according to their covariance matrices and repeating 
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TABLE I. Reconstruction efficiencies, including B{K° Kg vr^vr ), and signal yields for 
channels with r/c- 



Channel Efficiency (%) Fit yield (events) 

B+ ^ VcK+,7]c ^ 00 22.0 lAtl'l 

B° 7]cK°, 7]c K°KTr 3.9 7.5t|^ 

B° ^rj,K\7],^ c^qy 6.2 l.OlJ;^ 



the fit procedure. Systematic errors due to uncertainty in the reconstruction efficiency are 
quoted separately. An additional source of uncertainty is the rjc branching fractions; B{ric 
K°Kti) has a relative error of 30.9% and i3(?7c ~^ 00) has a relative error of 39.4%. Since the 
measurements of these branching fractions were made by experiments running at the J/ip 
mass, they have a common error of 28.3% due to the uncertainty of the branching fraction 
B{J/ip — >■ ^rjc). We quote this common error on the combined result as coming from the t]c 
branching fraction uncertainty. The remaining errors of 12.4% and 27.4% are included by 
smearing the likelihood functions for different sub-modes before they are combined. 

Assuming equal production of charged and neutral B mesons in T(4S') decay, the branch- 
ing fractions are B{B+ r]cK+) = (0.69l[]:i ± 0.08 ± 0.20) x lO'^ and B{B° t]cK°) = 
(1.091 o;42 ± 0.12 ± 0.31) X 10 ^, where the ffist error is statistical, the second error is from 
the uncertainty on the reconstruction efficiency and the third error is due to the uncertainty 
of B{J/il) 7'7c)- The statistical significances, with PDF shape uncertainties included, are 



5.2 standard deviations for the charged B decay and 4.8 for the neutral channel ||T2[. The 
confidence level of the ffis are 54% and 76% for the charged and neutral decay modes, respec- 
tively. The combined M{B) projections and the log-likelihood functions of the branching 
fractions are shown in Fig. [^. 

As a cross-check we have done a counting analysis using more stringent selection criteria 
\M{B) - 5.2801 < 0.007 GeV, \AE\ < 0.040 GeV and \M{r]^) - 2.9798| < 0.025 GeV. The 
results are statistically consistent with those stated above. We have checked the sensitivity 
of our result to the large spread of the rjc width measurements ||13| by repeating the analysis 
using values of 10 and 24 MeV for the rj^ width. The central values of the branching fractions 
are 0.64 x 10^^ and 0.78 x 10^^ for the charged decay and 1.02 x 10"^ and 1.30 x 10"^ for 
the neutral decay, correspondingly. 

The branching fraction for the decay B — > rjcK can be related to that of -B ^ J/ipK 
||T3| in the factorization approximation by taking into account the phase space difference 
and hadronic current dynamics The ratio of decay rates is proportional to the ratio 

of the decay constants squared: r(_B rjcK)/V{B — > J/ipK) = D^fr^J fj/^Y, where the 
coefficient D expresses the evaluation of the decay dynamics. We calculate the ratio of the 
decay constants from the measured branching fractions using weighted averages of charged 
and neutral modes and the theoretical estimates of D. Our results are consistent with the 
phenomenological expectations given in Table |I[ Using the J/ip decay constant evaluated 
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from dilepton rates, /j/^ = 405 ± 14 MeV and predictions of D by Ahmady and Mendel 
we obtain /^^ = 335 ±52 ±47 ±12 ±25 MeV, where the first error is due to the statistical 
and systematic errors on the exclusive branching fractions, the second error is due to the rjc 
branching fractions, the third error reflects the uncertainty in the J /ip decay constant, /j/,/,, 
and the last error is due to D. 



TABLE II. Theoretical estimates of T[B rj(.K) /T{B J/tpK). The first, second and 
fourth columns give phenomenological evaluations of different quantities. The third column lists 
our estimate of the decay constant ratio for a given model, where the first error originates from the 
branching fraction uncertainties and the second is due to the quoted error in D. 



D fyjfj/jj fyjfj/j, (exp.) r(B^J/W) 

^ 2.68 0.78 ±0.13 0.54 ± 0.10 ± 0.00 1.64 ±0.27 @ 
1.12 ±0.17 1.20 ±0.04 0.83 ±0.17 ±0.06 1.6 ±0.2 || 
1.11 ±0.15 0.99 0.84 ±0.17 ±0.06 [0.94,1.24] @ 
1.11 ±0.15 1.03 ±0.07 0.84 ±0.17 ±0.06 1.14 ±0.17 § 
1.43 ±0.29 0.81 ±0.05 0.73 ± 0.16 ± 0.09 0.94 ± 0.25 @ 



To search for the decay B — > XcoK we follow the same procedure as described above. Part 
of the background comes from B decays to J/ip or D mesons, such as 5 — > Dn, D Kit. In 
these cases the reconstructed B momentum is close to that of the signal and the background 
peaks in M{B). About one half of such background is removed by vetoing leptons and Kn 
combinations in the vicinity of charged and neutral D mesons masses. We do not model the 
behavior of the remaining background. Instead, we make a ±7 MeV cut on M{B). Only the 
AE and M{xco) variables are used to extract the signal. Mistaking pions from Xco decay as 
kaons shifts the B energy upwards by 120 MeV. The AE region is made asymmetric to keep 
event samples from different Xco decay channels from overlapping (—150 < AE < 60 MeV 
for the Xco K^K~ submode and —60 < AE < 150 MeV for the Xco tt+tt" submode). 
In spite of the lepton veto some of the background due to — > /z"*"//" decays remains in 
the sample because of restricted muon system acceptance. This background contributes 
to the M(xco) sideband of Xco K^K~ submode, hence we restrict the signal plane to 
-^(Xco) > 3.28 GeV. We remove any possible contribution from Xc2 tt+tt" or K^K~ by 
imposing a skew veto cut AE > M(xco) — 3.5 GeV. 

The final results were extracted from the limited region of AE and M(xco)- The same 
unbinned maximum likelihood procedure was used as for the B rj^K analysis. The 
observed yields are not statistically significant and the resulting 90% confidence level upper 
limits are 0.48 x 10"'^ and 0.50 x 10"'^ for the charged and neutral modes, respectively [p2 |. 

In summary, we have observed the decay B rjcK in both charged and neutral modes 
with branching fractions similar to those for B J/ipK. By comparing the rates of the 
decays with rjc and J/'?/', we have extracted the ?7c decay constant. The channel B° — >■ rjcK" 
can be used to extract the value of sin(2/3) via measurement of time-dependent asymmetry. 
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We have also set upper limits on S — > Xco^ decays that restricts a possible enhancement of 
the XcO production due to the color octet mechanism. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, the Natural Sciences and 
Engineering Research Council of Canada, the A. P. Sloan Foundation, the Swiss National Sci- 
ence Foundation, the Texas Advanced Research Program, and the Alexander von Humboldt 
Stiftung. 
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FIG. 1. Beam constrained mass spectra and the log- likelihood functions for charged (a and c) 
and neutral (b and d) B decay channels. The mass spectra include both rjc channels combined, with 
r]c — > ^-Iso shown separately as the shaded area. The solid line in figures (a) and (b) displays the 
signal plus background combined shape. The dashed line corresponds to the background shape only. 
A cut on the signal likelihood using all variables except M{B) is used to make these projections. 
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